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I ■ Abstract. 

, We study by density functional and large scale tight-binding transport calculations 

the electronic structure, magnetism and transport properties of the recently proposed 
graphene ribbons with edges rolled to form nanotubes. Edges with armchair nanotubes 
I present magnetic moments localized either in the tube or the ribbon and metallic or 

half-metallic character, depending on the symmetry of the junction. These properties 
. have potential for spin valve and spin filter devices with advantages over other proposed 

^ I systems. Edges with zigzag nanotubes are either metallic or semiconducting without 

O ■ affecting the intrinsic mobility of the ribbon. By varying the type and size of the 

nanotubes and ribbons offers the possibility to tailor the magnetic and transport 
properties, making these systems very promising for applications. 

^ I The atomic structure of graphene edges is important for the determination of 

I the electronic and magnetic properties of graphene, especially for narrow graphene 

2 : nanoribbons[Il |2l|3llll[5l[6l[71|8l|9l[ini[IIl[12]. Recent theoretical work[l3] on the 

' stability of different graphene edges structures has shown that graphene edges can fold 

^ ■ back on themselves and reconstruct as nanotubes, with low formation energy (see atomic 

structures in figured]). In this article, we show that, beside protecting the edges from 
^ . contamination and reconstructions, nanotubes at the edges may lead to magnetism and 

^ \ are not detrimental for the electronic mobility despite the row of sp^ hybridized atoms 

at the ribbon-tube junction. We study the electronic and magnetic properties of these 
systems by a combination of density functional theory (DFT) and large scale tight 
binding (TB) simulations of transport properties. Our calculations suggest that these 
systems could be used for a variety of applications that we sketch in figure |2j 

We consider systems formed by a nanoribbon terminated on both sides by the 
same armchair (AC) or zigzag (ZZ) nanotube. We notice that a ribbon with AC 
edges is terminated by ZZ nanotubes and a ribbon with ZZ edges is terminated by 
AC nanotubes. Nanoribbons terminated by AC nanotubes present interesting magnetic 
properties. By rolling the ZZ edges of a nanoribbon, two types of AC nanotubes can 
be formed, as shown in figure [H If the atoms at the nanoribbon ZZ edge scroll and 
bind to the same sublattice sites within the nanoribbon, the formed AC nanotube has 
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Figure 1. (a-h) Structure and spin density for symmetric (left) and asymmetric (right) 
AC nanotube terminated edges. (a,b) Atomic structure; (c,d) Side, (e,f) 3D and (g,h) 
top view of schematic spin representation. The box with dashed hne in (g,h) indicates 
half of the unit cell in the DFT calculation. For the symmetric case, the contribution 
to the magnetic moment from atoms belonging to the A- and B- sublattices are niA = 
0.958 ^j-B and itib = -0.217 (per half unit cell), and for the asymmetric case uia = 
1.048 fXB and tub = -0.300 /is. The maximum of the magnetic moment is located on 
the atom indicated by the cirrOW with the VclIuG ??lmaa; 

= 0.379 ^B in both symmetric 
and asymmetric cases. (i,j) 3D and top view of AC nanoribbons terminated by ZZ 
nanotubes. The box with dashed line in (j) indicates half of the unit cell in the DFT 
calculation. 
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Table 1. Total spin magnetization S (in unit oi fiB per unit cell) as a function of 
the index of AC nanotube {N, N) and size of nanoribbon P (in unit of ZZ rows). The 
results in this table are for samples obtained by rolling the edges of the same unrolled 
graphene ribbon with width of 40 ZZ rows. 
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mirror symmetry with respect to the nanoribbon plane; if the bonding sites belong to 
opposite sublattice, there is no such kind of symmetry (compare figure [T]d to figure 
[IK). We call these two cases symmetric and asymmetric which corresponds to armchair 
and armchair-like in Ref. [13], respectively. The common point of these two cases is 
that the sublattice symmetry is broken, because all the sp^ hybridized carbon atoms 
at the junction belong to one sublattice. Due to Lieb theorem [HI [12], this gives the 
possibility of spin polarization around the junctions. Since the theorem applies to the 
Hubbard model, accurate calculations for the real system are necessary to investigate 
this possibility. 

In order to study the magnetic properties, we performed spin polarized DFT 
calculations by SIESTA [TH [151 [iSj- We used generalized gradient approximation 
with Perdew-Burke-Ernzerhof parametrization (GGA-PBE) [T7| and a standard built-in 
double-^ polarized (DZP) [18] basis set to perform geometry relaxation. We found that, 
in both symmetric and asymmetric cases, there is spin polarization near the ribbon-tube 
junction, i.e., near the sp^ hybridized carbon atoms. Note that the bond distance of these 
four-fold coordinated atoms is 1.54 A like in diamond. The spin polarization is mainly 
located in the nanoribbon for the symmetric case and within the tube for the asymmetric 
case (see the isosurface plot of the spin density together with its symbolic representation 
in figure [It-h) . The up/down spins are distributed over the A/B sublattices respectively. 
The label nimax indicates the atom with the highest magnetic moment. 

For the symmetric case, the value of the spin polarization increases with increasing 
nanoribbon width, and saturates at 1.50 fis per unit cell, when the nanoribbon width 
is wider than 16 ZZ rows; For the asymmetric case, the spin polarization is always 1.50 
fiB per unit cell, irrespective of the nanoribbon width and nanotube radius (see Table 

For the asymmetric case the spins are located inside the two nanotubes and therefore 
the exchange interaction between opposite edges is negligible. For the symmetric case, 
the spins on the two edges are coupled antiferromagnetically, similarly to hydrogen 



Electronic, Magnetic and Transport Properties of Graphene Ribbons Terminated by NanotubesA 



a be 




Figure 2. Sketch of spintronics devices based on carbon nanoribbons terminated by 
AC nanotubes. (a) A spin valve based on the symmetric case: a gate can be used to 
switch from the antiferromagnetically coupled state to the ferromagnetically coupled 
excited state, favouring spin transport from one nanotube to the other across the 
ribbon, (b) For both symmetric and asymmetric case, high magnetorcsistance could 
be achieved by applying magnetic fields of different sign at the ends of the nanoribbons, 
as proposed in [19] for ZZ nanoribbons. (c) For the asymmetric case, a gate along the 
ribbon could be used to switch between two half-metallic energy regions to realize 
either a spin filter or a spin valve. 



terminated graphene edges [201 [2T| : for the structure shown in figure [T^, the energy of 
antiparallel spin configurations is 22 meV per unit cell (see figure [1^) lower than for 
parallel configurations. This sizeable coupling across the ribbon makes the symmetric 
systems promising as spin valve devices [22]. In figure |2^ we show a configuration 
similar to that proposed for dumbbell graphene structures on the basis of the Hubbard 
Hamiltonian in the mean field approximation |23]. A gate could be used to bring the 
system from the antiferromagnetically coupled state to the ferromagnetically coupled 
excited state, favouring spin transport from one nanotube to the other across the 
ribbon. Moreover, both for the symmetric and asymmetric case, the magnetic moments 
along the ribbon-tube junction are qualitatively similar to the case of ZZ edges of 
nanoribbons. Therefore, high magnetoresistance could be expected, as proposed in 
[To] for nanoribbons with ZZ edges, by applying magnetic fields of different sign at the 
ends of the nanoribbon. A sketch of this device for our systems is shown in figure [2b. 

The spin polarized density of states (DOS) reveals other features of interest for 
spintronics related to half-metallic character. In figure [3^ and [3b we show the spin 
polarized DOS for the symmetric and asymmetric case respectively. We see that the 
symmetric case is metallic for both spins in the whole range of energy. The asymmetric 
case, instead, is a half-metal near the Fermi energy Ep, namely it is metallic for spin 
up and insulating for spin down. The half-metallic character of our systems provides 
opportunities as spin filters without the need of external electric fields [21], magnetic 
fields [25]. ferromagnetic strips [26], impurities [27 [ [28 | [29] or defects [30] [3T]. Furthermore, 
there is the opposite half-metallic character at higher energies. Around 0.4 eV, there is 
insulating character for spin up and metallic character for spin down. As sketched in 
figure [2t, a gate along the ribbon could be used to switch between these two half- metallic 
energy regions and affect selectively the spin transport. 
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Figure 3. Density of states calculated by DFT for ZZ nanorribons with (a) symmetric 
and (b) asymmetric AC nanotubc terminated edges (A'' = 5 and P = 20 for both cases) . 
The arrows in (b) indicate the half-metalic energy regions for either spin up or spin 
down. 



We come now to the case with ZZ nanotubes shown in figure[T}-j. For this case, there 
is only one type of ribbon-tube junction that preserves sublattice symmetry implying 
that there is no magnetization nor midgap states [H [12]. The electronic structure and 
transport properties, however, strongly depend on the AC ribbon width and on the ZZ 
tube radius. In TB models, a AC nanoribbon is metallic if the number of AC rows 
is equal to 3/ + 2, where / is a positive integer, and semiconducting otherwise[Il 132] . 
Furthermore, ZZ nanotubes are metallic for index equal to 3/[33|. In more general 
models, the properties of AC nanoribbons and ZZ nanotubes may differ from the ones 
predicted by TB, due to possible self passivation of the edges for nanoribbons and for the 
cr — TT band mixing for small nanotubes [33]. By using DFT calculations, we found that 
our joined system becomes a semiconductor with a gap of the order of few hundreds meV 
if both nanoribbon and nanotube are semiconducting. The energy gap as a function 
of geometry is shown in figure HI The size of the nanotube has to be large enough for 
the opening of a band gap (figure H^.) . For the joined system with semiconducting ZZ 
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Figure 4. Band gap size-dependence in ZZ nanotube terminated AC nanoribbons. (a) 
Fixed nanoribbon with 15 AC rows for varying size of the terminating ZZ nanotube. 
(b) Fixed ZZ nanotube (8, and 7, 0) edges for varying width of the AC nanoribbon. 



nanotubes, there is a clear periodicity (3 ZZ rows) in the dependence of the energy gap 
on the nanoribbon width (figure Hb). For the studied cases, the value of the gap varies 
between 30 and 600 meV. 

Since we want to calculate the transport properties by means of a simpler model, 
suitable for large samples, we have also calculated the energy structure of our system by 
TT-band TB calculations where we consider only the nearest-neighbor hopping t = 2.7 
eV between the carbon atoms [3^ 135] for all bonds. The comparison between the 
gaps calculated by TB and DFT shown in figure HIi gives the same periodicity and 
a qualitative agreement. For the case of a AC nanoribbon with 15 rows terminated 
by (8,0) ZZ nanotubes we compare in figure [5^ the TB and DFT DOS, which again 
are in qualitative agreement as to support the validity of the transport calculations we 
show next. The electronic transport properties of a semiconducting nanoribbon with or 
without nanotube terminated edges are obtained by using large scale TB simulations 
with about thirty million carbon atoms [M| [35]. In figure [5]d we show that the electronic 
mobility parallel to the edges as a function of charge density are quite similar in these 
two cases. The mobility of the joined system is about 2250 cm?V~^s~^ at charge density 
Ue ~ lO^^cm^, which is only slightly smaller than the one of the AC nanoribbon at the 
same charge concentration. 

In summary, we have studied the electronic and magnetic properties of graphene 
nanoribbons with three types of nanotube terminated edges. The spin magnetization is 
found to be 1.5 /x^ per unit cell in the ground state of both symmetric and asymmetric 
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Figure 5. (a) Comparison of DOS of a nanoribbon with 15 AC rows terminated by 
(8,0) ZZ nanotubes calculated by DFT and TB. (b) Comparison of the mobility of a 
nanoribbon with 15 AC rows, terminated either by (8,0) ZZ nanotubes or by AC edges 
in TB. The charge density rig is obtained from the density of state p in panel (a) by 
Ue (E) = p (e) dt. The sample contains 600000 x 47 carbon atoms for the case with 
ZZ nanotube terminated edges, and 2000000 x 15 for the case with AC edges. 



AC nanotube terminated edges. For symmetric AC nanotube terminated edges, the spin 
density is located in the ribbon whereas, for the asymmetric case, it is located within the 
tube. In the ZZ nanotube terminated edges, there is a band gap opening of the order 
of few hundreds meV, if the constituent tube and nanoribbon are both semiconducting. 
The conductivity and mobility in the presence of ZZ nanotube terminated edges is 
comparable to the one of the AC nanoribbon itself. 

Our calculations suggest that these systems are not only advantageous because the 
edges are protected against any kind of chemically induced disorder but also because. 
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by tailoring the ribbon/tube structure, they offer a wealth of possible applications for 
transport and spintronics. 
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